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Abstract 

We consider the exterior Cauchy-Dirichlet problem for equivariant wave maps from 3 + 1 dimen- 
sional Minkowski spacetime into the three-sphere. Using mixed analytical and numerical methods 
we show that, for a given topological degree of the map, all solutions starting from smooth finite 
energy initial data converge to the unique static solution (harmonic map). The asymptotics of 
this relaxation process is described in detail. We hope that our model will provide an attractive 
mathematical setting for gaining insight into dissipation-by-dispersion phenomena, in particular 
the soliton resolution conjecture. 
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INTRODUCTION 



Dissipation of energy by dispersion is the key mechanism of relaxation to a static equi- 
hbrium in infinite dimensional Hamiltonian systems on unbounded domains It occurs 
in a variety of physical situations ranging from dynamics of gas bubbles in a compressible 
fluid [2] to formation of black holes in gravitational collapse jsj. Despite great physical 
importance, mathematical understanding of dissipation by dispersion phenomena is still in 
its infancy, especially in the non-perturbative regime. Thus, to make progress it would be 
useful to have a good toy-model of this phenomenon. Among the desired features of such a 
model we list the following: 1) the evolution equation is as simple as possible, preferably a 
scalar semilinear equation, 2) solutions are globally regular in time for all reasonable initial 
data, so that one need not to worry about the development of singularities, 3) the static 
equilibrium is unique (so that there is no competition between different attractors) and rigid 
(so that no modulation analysis is needed), 4) the linearization about the static solution has 
a purely continuous spectrum (no unstable, neutral or oscillatory modes), 5) the hnearized 
waves decay fast enough to start the Duhamel iteration argument (the main technical tool 
in proving asymptotic stability). If all these properties hold, then it is natural to expect 
that the static solution is a global attractor, i.e., all solutions converge to it as time goes 
to infinity while the excess energy is being radiated away to infinity. This is a particularly 
simple case of the soliton resolution conjecture which asserts that generic global in time 
solutions of nonlinear dispersive equations resolve asymptotically into a superposition of a 
coherent structure (solitons) and free radiation . 

The purpose of this paper is to propose a model which enjoys all the properties listed 
above. We hope that the insight gained by studying this model may help in understanding 
real physical problems. The rest of the paper is organized as follows. In section 2 we define 
the model. Static solutions are discussed in section 3. In section 4 we describe the spectral 
properties of the linearized problem. Finally, these analytic results are confronted with 
numerical simulations in section 5. 
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MODEL 



A map U : ^ M from a spacetime [M.,gai3) into a Riemannian manifold (A/", Gab) is 
said to be the wave map if it is a critical point of the action 



where x"" and X"^ are local coordinates on A4 and Af, respectively. Variation of the action 
dl]) gives a system of semilinear wave equations 

where := -^^d^ {g"^\/—g dp) is the wave operator associated with the metric gajs and 
are the Christoffel symbols of the target metric Gab- In this paper we consider the 
wave maps from M. = M^^^, the 3 + 1 dimensional Minkowski spacetime, into N' = S^, 
the unit 3-sphere. On M^+^ we use standard spherical coordinates {t, r, 6, 0) and on we 
choose polar coordinates X^ = (\E', 9, $), hence 

GABdX^dX^ = d^^ + sin^ ^ {dQ^ + sin^ 6 d*^) . (3) 

In addition, we assume that the map U is spherically £-equivariant, that is 

^ = uit,r), (e,0) = fi,(^,0), (4) 

where fli : S"^ ^ S*^ is a harmonic eigenmap with eigenvalue i{i + 1) (£ G N). Under these 
assumptions the wave map equation (|2]) reduces to the single scalar wave equation 

2 , + 1) 

— il + u" -\ — u — ^ — sin 2u = , (5) 

where dots and primes denote derivatives with respect to t and r, respectively. 

Eq.© is invariant under scaling: if u{t,r) is a solution, so is Ux(t,r) = u{t/\,r/\). The 
total conserved energy scales as E[ux] = XE[u], which means Eq.(|5]) is supercritical and 
singularities are expected to develop for some initial data. In fact, it is known that Eq.(|S]) 
has infinitely many self-similar solutions jsj which are explicit examples of singularities 
forming in finite time from regular initial data gJ. Moreover, the ground state self-similar 



solution is stable and acts as an attractor for generic large initial data 



Note also 



that the scaling symmetry precludes existence of nontrivial static solutions. The formation 



of singularities for generic initial data and the absence of nontrivial static solutions make 
Eq.(|5]) uninteresting from a physical point of view. One way to cure these deficiencies is to 
modify the action ([T]) by adding a suitable higher order term with subcritical scaling. Such a 
term breaks the scaling symmetry and prevents shrinking of solutions to zero size. The best 
known modification of this type is the Skyrme model [9;]- This model admits stable static 
solutions (so called skyrmions) and enjoys global regularity in time [lo|. Unfortunately, 
these nice features are achieved at the expense of a significant increase of complexity: the 
evolution equation is quasilinear and very difficult to analyze (note that it took fifty years to 
prove that no singularities can form during the evolution). The numerical evidence for the 

n 

soliton resolution conjecture in the Skyrme model was given in but to our knowledge 
there is no corresponding rigorous result. 

Motivated by a desire to have a simpler model than Skyrme's and yet possessing similar 
properties we proceed in a different manner. We keep the action unchanged but modify the 
domain of the wave map. Namely, we consider Eq.(l5]) outside a ball of radius i?, that is for 
r > R and impose the Dirichlet condition at the boundary u{t,R) = [12]. Hereafter, we 
choose the unit of length such that R = 1. We are interested in the long-time behavior of 
solutions starting from smooth initial data with finite energy 

1 1"^ 

E{u) := - / (r^M^ + r\''^ + i{i + 1) sm\) dr < oo . (6) 

Note that finite energy solutions must satisfy the boundary condition at infinity u{t, oo) = 
A^vr, where an integer N is the topological degree of the map (sometimes called the winding 
number). Since the degree cannot change during the evolution, the initial value problem 
for Eq.(j5]) breaks into infinitely many disjoint topological sectors labeled by the degree N. 
Note also that global regularity in time is automatically guaranteed because by spherical 
symmetry singularities for Eq.([5]) cannot occur for r > 0. In other words, by cutting off the 
ball from the domain of the map we make the equation effectively subcritical. 

Knowing that solutions exist for all times, it is natural to ask how they behave as t — t- oo. 
Below we show that for each £ G N and each G Z there exist a unique static solution 
ui^Nij) which is linearly stable. We conjecture that the solution U(,^n is a universal global 
attractor, that is all solutions starting with smooth initial data of degree N asymptotically 
converge to it. The numerical evidence supporting this conjecture is presented in section 5. 
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STATIC SOLUTIONS 



To study static solutions it is convenient to use new variables x = Inr and v{x) = u{r). 
Then, time-independent Eq.dS]) takes the autonomous form 

v" + v' - ^ ^ ' sin 2v = 0, (7) 

which in mechanics describes the motion of a pendulum with constant friction. A regular 
finite energy solution of degree N corresponds to a trajectory which in the phase plane 
{v, v') starts from a point (0, a) and goes to the point (A^vr, 0) as x — )■ oo. The existence and 
uniqueness of such a trajectory for each N follows from elementary phase-plane analysis. 
The phase portrait displaying the flow of the vector field X = {v', —v' + + 1) sin 2v) 
is depicted in Fig. 1. The flow has stable focal points at ((A^ + l/2)7r,0)) with eigenvalues 
a = (—1 ± i ^/A£{£^^l)'^^) / 2 and saddle points at (A^tt, 0) with eigenvalues a = i and 
a = —i — 1. From each saddle point there emerge two separatrices corresponding to the 
stable [a = —i — 1) and unstable {a = €) manifolds. It is routine to show that close to the 
saddle the stable manifold can be parametrized in the form 

v{x) = NT:-h e-^'+^> + O (e-3(^+i)-) , (8) 

where the coefficients of all higher order terms are determined by the parameter h. Gener- 
ically the point (0, a) belongs to the basin of attraction of one of the focal points but for 
isolated values a = ag^N this point lies on the separatrix which runs into the saddle {Nn, 0) 
along the stable manifold. 

Translating the above analysis back to the original variables we conclude that for each 
integer N there exists a unique static solution ug^Nij") of Eq.Q satisfying the boundary 
conditions M^^iv(l) = and u^7v(oo) = Nti. This solution is parametrized by the pair of 
numbers (a£,Ar, &^,Ar) such that 'u^^(l) = a^ ^r and 

M,,^(r) = ATtt - ^ + 0(r-=^('^+i)) as r^oo. (9) 

The parameters a^^N and hi^^ can be determined numerically with the help of a shooting 
method. A few values of these parameters for different i and N are listed in Table I, while 
Fig. 2 depicts several profiles ui^iy{r). Let us note that the static solutions described above 
were first found by Balakrishna et al. pJi] who used them as a simple model of nucleons. 
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FIG. 1: The phase portrait for Eq.Q for £ = 2. 



1 = 1 1 = 2 1 = 2, £ = 4 i = 5 1 = 6 



be,i 


3.786299 4.327397 4.800690 5.226581 5.617059 5.979801 
4.847841 6.147165 7.240761 8.203436 9.073200 9.872569 


ae,2 
be,2 


6.983263 7.601218 8.161078 8.677933 9.161071 9.616701 
15.810784 29.162429 46.410705 67.463775 92.172786 120.372175 



TABLE I: Parameters of the static solutions. 



Remark 1. It does not seem feasible to express the solutions U£,Ar(r) in closed form, however 
for large i one can obtain a good analytic approximation for ue^i as follows. Let y = {i + l)x 
and v{y) = v{x). Then, in the limit £ — )■ oo Eq.(l7]) becomes v" — | sin 2v = 0. The solution of 
this limiting equation satisfying the boundary condition v{oo) = tt is v{y) = 2 arctan (Ae^), 
where A is an integration constant. Returning to the variable x and setting A = 
(to agree with ([8])), we obtain the large i approximation v{x) ~ 2 axctan (^■^^e^^'^^''^^ , 
which is valid outside the boundary layer of width 0{i~^) near x = 0. Using matched 
asymptotics one can systematically match this outer solution to the inner solution satisfying 
the boundary condition v{0) = 0,f(0) = a^^i, however for practical purposes the following 
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ad hoc modification (which does not alter the large x asymptotics) is good enough 

veAx) ^ ve,i{x) = 2 arctan A(e(^+i)- - e-('+')")^ . (10) 
As shown in Fig. 2, this approximation is quite accurate already for i = 1. 




FIG. 2: Left: Profiles of the static solutions for iV = 1,2 and £ = 1,2,3. Right: Deviation of the 
approximate solution ue^i{r) = V£^i{lnr) from the true solution ue^i{r). 

LINEARIZATION ABOUT STATIC SOLUTIONS 

Substituting u{t, r) = ui^]\r{r) +w{t, r) into Eq.(j5]) and linearizing we obtain the evolution 
equation for small perturbations around the static solution 

iij-^'- Inj' + ^^^^ w + V,^j,{t)w = , V,Mr) = -^^il±^ sin^ ue^r) ■ (H) 

Assuming harmonic time dependence, w{t,r) = e~*^*^/'(r), we get the eigenvalue problem 

- ^" - l^' + ^^7^ ^ + ^^Ar)^ = (12) 

for ip E L2 ([1, 00), r^dr) and satisfying ip^l) = 0. The spectrum is purely continuous > 0. 
To see this, note that ^{r) = r^^^^^^ satisfies Eq.(^ for = 0, as can be readily checked 
by differentating the time-independent Eq.(l5]) with respect to r. The presence of this zero 
mode is the consequence of scale invariance of Eq.(|5]). The zero mode is not an eigenfunction 
because it does not vanish at r = 1, nonetheless the fact that it has no zeros (since M£,Ar(r) is 
monotone) implies by the Sturm oscillation theorem that there are no negative eigenvalues. 
Therefore, the static solutions u^^n are linearly stable. 
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FIG. 3: The potentials Vi^N{r) for N = 1,2 and £ = 1,2, 3. 



The hnear stabihty is only the first necessary step in estabhshing asymptotic stabihty. 
To understand the latter we need to determine the quasinormal modes (also called complex 
resonances or scattering frequencies), that is solutions of Eq. fll2p which vanish at r = 1 and 
satisfy the outgoing wave condition tp{r) ~ ^e*^** with Q'(A) < for r oo. The fundamental 
(that is, the least damped) quasinormal mode is expected to govern an intermediate phase 
of the relaxation to the static solution. 

Let us first describe the quasinormal modes about the trivial solution ue^ = 0. In this 
case Eq.f ll2p takes the form of the spherical Bessel equation 

- r - 1^' + ^^^^ = , (13) 

whose exact outgoing solutions are given by the spherical Hankel functions of the first kind 
iplr) = h^p{\r). Thus, the quantization condition for the eigenfrequencies of quasinormal 
modes is '?/'(l) = h^p{X) = which is a polynomial equation of order hence has exactly £ 
roots. For example, using the formulae for the first two spherical Hankel functions 

yiz-^ z J yiz'^ z^ iz J 

one gets A = — i for ^ = 1 and A = ±-\/3/2 — 3i/2 for £ = 2. From the known results on 
zeros of Hankel's functions it follows that the eigenfrequency of the fundamental quasinormal 
mode for large (. is (see Proposition E.2 in {2!) 

X, = u-Cu-^+ ^C,^y-^ + O , (15) 

where z/ = £ + 1/2, C = 2~^e~^r]i, and r/i ~ —2.338 is the first zero of the Airy function. 
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For 7^ it is not possible to calculate the quasinormal modes analytically, e.g. because 
the potential Vg,Ar(r) in f|T2|) is not known in closed form, thus we must resort to numer- 
ical methods. We use a shooting method which goes as follows. First, we introduce new 
independent and dependent variables 

p = - , 7]{p) = re-'^'ipir) . (16) 
r 

The advantage of this transformation is two-fold: i) it peels off the oscillatory behavior of 
ilj{r) at infinity and ii) compactifies the semi- infinite interval 1 < r < oo to the finite interval 
< p < 1 (switching the endpoints). In terms of the new variables Eq.( [T2|) becomes 

p^r^" + (2p - 2iX)r]' - i{£ + l)r/ + 2i{£ + 1) sin^ (m^^^) t] = . (17) 

Next, we solve this equation numerically away from the endpoints starting from the desired 
"initial values". At p = 1 the initial value is ri{l) = (and an arbitrary nonzero ^^'(1)). At 
p = the desired outgoing solution is defined by the asymptotic expansion 

£ ■ oo 

"W-Ep||^^> E '.V. (18) 

The first finite sum in ( ITSl) corresponds to the free outgoing solution. The second infinite 
sum is the divergent asymptotic series (since p = is the irregular singular point) whose 
coefficients Cj^X^h^jq) can be calculated recursively using Q. We use the series (ITSl) to 
impose the initial value at p = 0.1. Note that the unwanted ingoing solution is of the order 

2|a(A)| 

0{e p ), hence it is difficult to keep track of it numerically, especially if |Q'(A)| is large. 
To suppress the contamination of the initial value at p = 0.1 by this exponentially small 
ingoing admixture it is crucial to compute the series (JTSll with very high precision. We 
achieve this precision by the optimal truncation of the asymptotic series (alternatively, the 
Fade summation could be used). The requirement that the logarithmic derivatives of the 
two numerical solutions launched from p = 0.1 and p = 1 join smoothly at a midpoint, 
say p = 0.5, serves as the quantization condition for quasinormal modes. The fundamental 
eigenfrequencies determined in this manner are listed in Table II. We denote the j-th (j = 
0, 1, . . . ) quasinormal eigenfrequency and eigenmode by = fi;^^ — iF;^^ and iplj^{r). 
We adopt the convention that for given A^ and I the damping rate F;^^ increases with j, 
hence the j = mode is fundamental. Below, to avoid notational clutter, we drop the index 
j = on the fundamental quasinormal modes. 
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Note that as ^ increases, the frequencies fi^^A? increase while the damping rates F^^at 
rapidly decrease to zero. This behavior was intuitively to be expected in view of the fact 
that the potential Vf,Ar(r) becomes narrower and deeper as £ grows (see Fig. 3) (hence the 
waves get increasingly squeezed and trapped as ^ grows). 





e = i 


£ = 2 


£ = 3 


£ = 4 


£ = b 


£ = % 






0.425517 


0.910515 


1.346222 


1.752249 


2.129364 


2.475564 




,1 


0.347445 


0.246627 


0.151154 


0.080210 


0.035763 


0.012945 




?,2 


0.228006 


0.436739 


0.632007 


0.814959 


0.982609 


1.137494 


r^ 


,2 


0.121276 


0.058779 


0.021354 


5.2862 X IQ-^ 


8.6491 X 10~^ 


1.0037 X 10"^ 



TABLE II: The fundamental eigenfrequencies A^^^r = ^t,N — i^£,N for = 1, 2. 



1 

10-2 
10-^ 
10"^ 
10-8 
10-10 

-6 -4 -2 2 4 6 

FIG. 4: The distribution of fundamental eigenfrequencies Xi^jy for A = 1, 2 and £ = 1, 2, . . . , 15 in 
the complex A-plane. Numerics indicates that n ~ £1/2 and P/ ^ ~ e ^ ^ > foi large £. 

To reassure ourselves and verify independently the above computation we solved nu- 
merically the linearized evolution equation (ITT]) . According to the scattering theory, for 
intermediate times the solution at a fixed point is space can be represented by the sum of 
quasinormal modes 

w{t, r) -^Yl Aj exp (-iA^^^t) ipl^ir) , (19) 
j 
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where the coefficients Aj depend on initial data. After a short transient time, all the quasi- 
normal modes but the fundamental one die out (see Fig. 5), hence we can readily determine 
the fundamental eigenfrequency by fitting the exponentially damped sinusoid to the numer- 
ical solution. The results of such a fit agree to four (for i = 1,2) and six (for ^ > 3) decimal 
places with the values given in Table II. To compute the parameters of excited modes we 

n 

applied the Prony method [13]. The outcome of this computation suggests that there are 
exactly quasinormal modes about the solution ue^N- 




20 40 60 80 

t 

FIG. 5: Quasinormal ringdown to zero for the solution of the linearized equation (llip for N = 2 
and £ = 10. After the passage of the direct signal from the initial data, an observer sitting at r = 5 
sees two clearly pronounced trains of exponentially damped oscillations which correspond to the 
first excited quasinormal mode with A}g 2 = 6.2925 — i 0.1718 and the fundamental quasinormal 
mode with Aio,2 = 1.6862 — z 9.1606 x 10^^°, respectively. 



Remark 2. For very late times the linear evolution described by Eq. flTT]) changes character 
from the quasinormal ringdown (which becomes negligible) to the polynomial decay called 
a tail. The tail is due to the backscattering of waves at large distances and its decay rate 
is insensitive to the presence of the ball. According to the linear scattering theory (see e.g. 
jl^). for the potential falling off as (/3 > 3) for large r and compactly supported initial 
data, the tail decays as t~'^, where 7 = 2£ + /3. In the case at hand, it follows from and 
(11 ip that V£^Ar(r) ~ ^-(2^+4) ^ hence 7 = 4£ + 4. We note for completeness that 

for the initial data which are not compactly supported and fall off polynomially at spatial 
infinity there is also a direct polynomially decaying in time signal from the data. 
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NUMERICAL EVIDENCE 



In this section we solve numerically Eq.(|5]) with the Dirichlet boundary condition u{t, 1) = 
and smooth finite energy initial data of various topological degrees N. The numerical 
technique is standard. We use the method of lines with a fourth-order Runge-Kutta time 
integration and fourth-order spatial finite differences. Our computational domain is large 
enough so that its outer boundary is causally disconnected from the region of interest. To 
ensure numerical stability, the high frequency modes were damped using the Kreiss-Oliger 
dissipation. To suppress round-off errors, all simulations were performed in the quadruple 
(128-bit) arithmetic precision. It is worth stressing that the Dirichlet boundary condition 
and smoothness imply that the Taylor series expansion at r = 1 takes the form 



u 



(t, r) = a{t){r - 1) - a{t){r -1)^ + ((r - l)^) 



(20) 



where a free function a{t) determines recursively all the coefficients of the series. Taken 
at t = 0, the expansion (120|) and its time derivative express the compatibility conditions 
between the initial values and the Dirichlet boundary condition. 



r 
10 



100 



27r 




H 1 — I I I I I I I 1 1 — I I I I I I 

t=8 - 



H 1 — I I I I I I I 



H 1 — I I I I I I 

t=16 - 



H 1 — I I I I I I I 



H 1 — I I I I I I 

t=32 - 




FIG. 6: Convergence to the static solution ni^i(r) (dashed line). 



Our main observation (illustrated in Fig. 6) is that for all i every solution of Eq. sat- 
isfying u{t,l) = and starting from smooth finite energy initial data of degree N converges 
(in any compact interval) to the static solution ue^N as t — > oo. A rigorous proof (and a 
precise formulation) of this claim is left as a challenge to the PDE experts. 
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Below we present some more quantitative aspects of the relaxation process. We focus on 
pointwise convergence in the timelike direction, that is we look at the time development of 
u{t, r) at a fixed point in space. We consider two kinds of initial data: (i) spatially localized 
perturbations of static solutions and (ii) arbitrary initial data with different next-to-leading 
order asymptotic behavior at spatial infinity than static solutions. We illustrate the case (i) 
with initial data of the form of the "kicked" static solution 

(r - ro)^~ 



m(0, r) 



2a2 



(21) 



with parameters e = 0.8, vq = 1.5, a = 0.1, and the case (ii) with initial data 

2fr- 1)2 



u{0, r) 
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arctan(r — 1), ^(0, r) = 



(22) 



1 + (r 

In both cases the intermediate phases of evolution are similar and have the form of quasi- 
normal ringdown, which is dominated by the fundamental quasinormal mode. Thus, this 
phase is well described by the linearized approximation (see Figs. 7,8,9). 
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FIG. 7: Pointwise convergence to the static solution for i = 1,2,3 and = (left panel) and 
= 1 (right panel) for exponentially localized perturbations of the static solution of the form 
(j2ip . Quasinormal ringdown for intermediate times yields to a polynomial tail for late times. In 
both cases the tail decays as (for £ = 1) and as (for i > 2). Note that for > 1, as 

I increases, the duration of the ringdown phase rapidly increases because the damping rate of the 
fundamental quasinormal mode decreases while the decay rate of the tail increases. 



The late time behavior is more complicated because there are three competing contri- 
butions to the tail: two kinds of linear tails mentioned in Remark 2 and the nonlinear tail 
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FIG. 8: Snapshots from the evolution starting from initial data ()2ip for i = 3, N = 1. Behind 
the wavefront, the spatial profile of the solution approaches the spatial profile of the fundamental 
quasinormal eigenfunction A 3? (e~*'*'^'^^*~*°^^/^3^i (r)) . The parameters A and tQ were fitted at t = 4.5 
once for all times. For better visibility the profiles are multiplied by e^^-^* (upper plots) and 
re^s,i{t-r) (bottom plot). 



caused by self-interaction of the field. The latter was described in detail in [l5|] in the frame- 
work of weakly nonlinear perturbation analysis, in particular it was shown there that the 
nonlinear tail decays as for £ = 1 and as for i > 2. For the initial data (121 p 

this tail dominates the late time evolution (Fig. 7) while for the initial data (1221) the linear 
tail coming from the deviation of asymptotic behavior of the data from the static solution 
is dominating (Fig. 9). 
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FIG. 9: Convergence to the static solution U£^i\f{r) for I = 1, 2, 3, 4 and = 1 for initial data (j22p . 
The deviation of these data from the static solution is of the order 0{r~^^^^^) for large r which 
leads to the linear tail decaying as t~'^ with 7 = 2 = 1), 7 = 4 = 2), 7 = 4 = 3), and 7 = 6 
= 4). The nolinear tail is subdominant. 

Acknowledgments 



We thank Patryk Mach for helpful remarks. This work was supported in part by the 
NCN grant NN202 030740. 



[1] A. Soffer, Dissipation Through Dispersion, in CRM Lecture Notes, Nonlinear Dynamics and 

Renornmalization Group, I.M. Sigal, C. Sulem (eds.). Vol. 27, (2001) 
[2] A.M. Shapiro, M.I. Weinstein, Radiative Decay of Bubble Oscillations in a Compressible Fluid, 

SIAM J. Math. Anal. 43, 828-876 (2011) 
[3] K.D. Kokkotas, B.C. Schmidt, Quasi-Normal Modes of Stars and Black Holes,. Living Rev. 

Relativity 2, 2 (1999) 

[4] T. Tao, Global behaviour of nonlinear dispersive and wave equations. Current Developments 

in Mathematics 2006, International Press. 255-340. 
[5] P. Bizoh, Equivariant Self-Similar Wave Maps from Minkowski Spacetime into 3-Sphere, Com- 

mun. Math. Phys. 215, 45 56 (2000) 
[6] J. Shatah, Weak solutions and development of singularities of the SU(2) a-model. Comm. 



15 



Pure Appl. Math. 41, 459469 (1988) 
[7] P. Bizon, T. Chmaj, Z. Tabor, Dispersion and collapse of wave maps, Nonlinearity 13, 
14111423 (2000) 

[8] R. Donninger, On stable self-similar blow up for equivariant wave maps, Comm. Pure Appl. 

Math. 64, 1029-1164 (2011) 
[9] T.H.R. Skyrme, A Non-Linear Field Theory, Proc. R. Soc. A 260, 127 (1961) 
[10] D. Li, Global wellposedness of hedgehog solutions for the (3+1) Skyrme model, preprint (2011) 
[11] P. Bizon, T. Chmaj, A. Rostworowski, Asymptotic stability of the Skyrmion, Phys. Rev. D 75, 
121702 (2007) 

[12] Searching the Uterature we found that this modification was first proposed by B. S. Balakr- 
ishna, V. Sanyuk, J. Schechter, and A. Subbaraman, Phys. Rev. D 45, 344 (1992) (as a model 
of the short-distance cutoff in QCD) but we are not aware of any further studies of this model. 

[13] G. Majda, W. Strauss, M.S. Wei, Computation of exponentials in transient data, IEEE Trans. 
Antennas and Propagation 37, 1284 (1989) 

[14] E.S.C. Ching et al.. Wave propagation in gravitational systems: Late time behavior, Phys. Rev. 
D 52, 2118 (1995) 

[15] P. Bizon. T. Chmaj, A. Rostworowski, S. Zajac, Late-time tails of wave maps coupled to 
gravity, Class. Quantum Grav. 26, 225015 (2009) 



16 



